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Why is a noble metal catalytically active? The role of the O-Ag interaction in the 

function of silver as an oxidation catalyst 
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Extensive density-functional theory calculations, and taking into account temperature and pres- 
sure, affords a comprehensive picture of the behavior and interaction of oxygen and Ag(lll), and 
provides valuable insight into the function of silver as an oxidation catalyst. The obtained phase- 
diagram reveals the most stable species present in a given environment and thus identifies (and 
excludes) possibly active oxygen species. In particular, for the conditions of ethylene epoxidation, 
a thin oxide-like structure is most stable, suggesting that such atomic O species are actuating the 
catalysis, in contrast to hitherto proposed molecular-like species. Copyright (2003) by The America! 
Physical Society. 
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The importance of reactant bond strength in hetero- 
geneous catalysis is demonstrated by the so-called "vol- 
cano curves"i2i^ which relate the reaction rate to the 
adsorption energy of the reaction intermediates. Typi- 
cally a maximum occurs in the rate at moderate values. 
This behavior underlies the expectation that a good cat- 
alyst should readily dissociate adparticles but not bind 
the fragments too strongly - a concept put forward by 
Sabatier in the early 1900's, which recent theoretical and 
experimental studies have re-addressed. 4,5 Despite this 
expectation, silver, as a noble metal, which binds ad- 
particles only relatively weakly on the surface, is a very 
important heterogeneous catalyst for various oxidation 
reactions. For example, partial oxidation of methanol 
to formaldehyde, as carried out at atmospheric pressure 
and temperatures of 800-900 KjSil as well as the selec- 
tive oxidation of ethylene to epoxide, also conducted at 
atmospheric pressure and at temperatures of 500-600 K. 8 

Although numerous investigations aimed at elucidat- 
ing the active oxygen species involved in the above- 
mentioned reactions have been carried out, it can be 
said that presently there is considerable confusion, and 
no unambiguous identifications have been made. For 
example, unclear issues concern the debate of molecu- 
lar (ozone-like )» versus atomioifi oxygen as the active 
species for ethylene epoxidation, the proposed impor- 
tance of "bulk-dissolved" oxygen and a "strongly bound, 
surface-embedded" species (which desorbs at 900 K) for 
the partial oxidation of methanol^ as well as the pos- 
sible role of chemisorbed or surface oxide-like species. 
Part of the reason for this lack of understanding is due to 
the difficulties associated with the pressure and materials 
gaps which exist for the oxygen-silver system. That is, 
the behavior is different under ultra-high vacuum (UHV) 
conditions, where the system can be analysed quantita- 
tively, to that under the high temperature and pressure 
conditions of catalysis, where it is much more difficult to 
obtain the same level of microscopic information. 

Our theoretical strategy, aimed at gaining an under- 



standing of the function of silver as an oxidation catalyst, 
in spite of its relatively weak binding capability, and to 
identify (and exclude) possibly active species, is to carry 
out extensive and systematic density functional theory 
(DFT) calculations for all conceivably relevant oxygen 
species. We take the effect of temperature (T) and pres- 
sure (p) into account through the oxygen chemical poten- 
tial, and thus determine the (T, p) phase-diagram which 
describes the surface phases from UHV right up to the 
conditions of real catalysis. We study oxygen at the (111) 
surface of silver as experiments indicate that this orienta- 
tion is an important crystal face for real silver catalysts 
in that at high temperatures, facets with this face re- 
sultAii Our general findings, however, are expected to 
be relevant for silver per se, and possibly to have impli- 
cations for gold, also a noble metal oxidation catalyst 
which, for the (111) surface, exhibits a restructuring and 
chemisorption of oxygen atoms at elevated temperatures 
(500-800 K) and atmospheric pressures; 13 i 14 similar to 
Ag(lll). They could also have consequences for copper, 
which catalyses the oxidation of methanol to formalde- 
hyde 

Our calculations employ the pseudopotentiaiiSiii 
plane-wave method^ with the generalized gradient ap- 
proximation (GGA)iS for the exchange-correlation func- 
tional. Five Ag layers are used in the supercell to model 
the O/Ag surfaces, which are created on one side of the 
slab, with a vacuum region of 15 A . The same k-point 
sampling is used for all structures. It corresponds to 21 
points in the irreducible part of the Brillouin zone of the 
(1 x 1) Ag(lll) surface cell. An energy cutoff of 50 Ry 
is used and full relaxation of the top two or three silver 
layers, and the O atoms, is taken into account. We in- 
clude the spin-polarization energy in calculation of the 
total energy of the free atoms and molecules. Further 
details are reported in Ref. [2(]. 

The oxygen species investigated, for a wide range of 
coverages, include: on-surface sites, surface- and bulk- 
substitutional sites, interstitial sites under the surface Ag 




FIG. 1: Selected results of average adsorption energies (with 
respect to atomic O) versus coverage for various structures: 
pure on-surface oxygen in fcc-hollow sites (red diamonds), 
surface- and bulk-substitutional adsorption (blue triangles), 
pure sub-surface oxygen in octahedral sites under the first 
("subsurface-octa.") and second ( "bulk-octa." ) Ag layers 
(black circles), an 03-like molecule adsorbed at a surface Ag 
vacancy (left-pointing black triangles), and oxide-like struc- 
tures (red squares). Oxygen atoms are depicted as yellow 
spheres and Ag atoms as the larger blue (and grey) spheres. 

layer and deeper in the bulk, as well as oxide-like struc- 
tures, and a molecular ozone-like species adsorbed at a 
surface Ag vacancy. 

The main important adsorption energies are summa- 
rized in Fig. 1, where for each "type" of oxygen atom, 
only the lowest energy structure is shown. For on-surface 
adsorption, the fcc-hollow site is preferred for all cover- 
ages investigated. For sub-surface adsorption, the octa. 
site under the first Ag layer is the most favorable. This 
is also the case for adsorption under the second Ag layer 
which may be taken to represent "bulk dissolved" oxy- 
gen. For all sub-surface sites, adsorption under the sec- 
ond layer is less favorable than under the first layer. In 
view of the unfavorable energy of oxygen in the bulk, high 
concentrations of bulk-dissolved oxygen in the perfect fee 
lattice of silver are improbable. 

Surface- or bulk-substitutional adsorption is seen from 
Fig. 1 to be unfavorable, as is the ozone-like species. To 
investigate the formation of oxide-like structures, we per- 
formed systematic investigations for geometries involving 
coverages of 0.50, 0.75, 1.0, and 1.25 ML using (2 x 2) 
cells. In particular, we tested combinations of on-surface 
fee and hep sites, and the three possible sub-surface sites 
under the first Ag layer. The atomic geometry of the fa- 
vored structure consists of 0.25 ML of oxygen on the sur- 
face and 0.5 ML in between the Ag layers. For increasing 
O concentrations, this structure "grows" by filling deeper 
lying sites between the Ag layers. The structure for a to- 
tal coverage of 1.25 ML oxygen is sketched in the inset 
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FIG. 2: Surface free energies for various low energy struc- 
tures as a function of the O chemical potential where A/io = 
MO - |£o 2 tal - For pressures of 10" 12 atm. (UHV) and 1 atm. 
(condition of catalysis), the corresponding temperatures are 
given. The temperature labels correspond to the tick marks 
and labels for the chemical potential. The labels "0.75 ML", 
"1.25 ML", and "2.25 ML", indicate the O-concentrations 
in the corresponding oxide-like structures. Accordingly the 
(4 x 4) structure has the label 0.375 ML, and the adsorption 
structures 1/4, 1/9, and 1/16 ML. At the top of the figure, 
the "material type" which is stable in the corresponding range 
of chemical potential are listed and indicated by the colored 
regions. 

of Fig. 1. In addition, we carried out calculations for the 
recently proposed (4 x 4) thin "surface-oxide" structure^ 
which involves an O-Ag-0 bonding unit (and electronic 
structure) similar to the upper tri-layer of the oxide struc- 
tures described above (see insets of Fig. 1). The average 
adsorption energies of all these oxide-like structures are 
displayed in Fig. 1, shown as red squares. Their atomic 
geometries are in fact very similar to the (111) surface of 
Ag 2 0. 

We turn now to investigate the effect of pressure and 
temperature on the stability of the various structures. To 
do this we calculate the surface free energy, 

~/(T,p) = (G-N Ag » As -No»o)/A , (1) 

where G = Gq^^^ — and the first and sec- 

ond terms on the right-hand-side are the free energies of 
the O/Ag surface under consideration and that of the 
reference system, i.e. the clean Ag(lll) slab, respec- 
tively. A is surface area in A 2 and A*o is the number of 
oxygen atoms. N Ag is the number of silver atoms (with 
respect to the reference system) and /ZA g is the Ag chem- 
ical potential, which is the free energy of an Ag atom 




in bulk fee silver. (Consideration of N& g and p,Ag are 
necessary only for the substitutional structures and the 
(4 x 4) phase) . The T and p dependence is given by , 
the oxygen chemical potential*" 

Ho<?,p) = l/2[^° 2 tal + Mo 2 (T,p°) + k B T In (^)],(2) 

where p° corresponds to atmospheric pressure and 
p,O 2 (T,p ) includes the contribution from rotations and 
vibrations of the molecule, as well as the ideal-gas en- 
tropy at 1 atmosphere. Here we use experimental values 
from thermodynamic tables. 

Through consideration of the difference in the vi- 
brational and entropic contributions of Go/Ag(m) an< ^ 
C?Ag(iii)' we find that it can be neglectediS? Thus, for 
evaluating the difference of the two G slab values, we re- 
place G slab by the total energy. We choose the zero of the 
O chemical potential to be half the total energy of O2 cal- 
culated as l/2££f 1 - E%£ Q - 2E^ - H f Ag?0 _ hulk . 
That is, as the energy difference between bulk silver ox- 
ide (^a^o) and bulk silver (2£^ al ), and the experimen- 
tal value of the heat of formation (-ffj^o-buik) which is 
-0.323 eVM 

The obtained results are displayed in Fig. 2. The 
range of A^o (defined as — ^^oT" 1 ) between the ver- 
tical dashed lines corresponds to the heat of formation 
of Ag20, i.e., to the range in which bulk silver oxide is 
stable. It can be seen that for more O-rich conditions 
(right) the thicker oxide-like structures are favored over 
chemisorption on the surface. Clearly, in the limit of 
an infinitely thick film, the corresponding line would be 
very close to the vertical one at A/io = —0.323 eV. It can 
be obtained from Fig. 2 that bulk silver oxide is stable 
to only around 350 K at atmospheric pressure, which is 
less than the experimental value of 460 K. At 460 K, the 
bulk oxide is unstable when the pressure is lower than 
10 3 atm. i.e. compared to 1 atm. as observed experi- 
mentally. This difference may be due to systematic errors 
of the DFT approach^ and/or to neglect of the entropy 
contributions. Therefore in Fig. 2 (and Fig. 3), we gener- 
ally expect that the temperature is underestimated, but 
the error is less than 110 K, while correspondingly, the 
pressure is overestimated. For the temperature values 
quoted below, we give an estimate of the upper limit in 
brackets, that is, 110 K higher than that obtained from 
Fig. 2. In spite of these uncertainties, we believe that the 
relative stability of the various systems and our general 
understanding and conclusions will not be affected. 

Due to the above-mentioned low thermal stability of 
silver oxide and the thicker oxide-like structures, our re- 
sults can safely rule them out as playing an important 
role in the catalytic reactions. For values of the O chem- 
ical potential further to the left in Fig. 2, and atmospheric 
pressure, the (4 x 4) phase is the most stable for the tem- 
perature range 350-530 K (incl. corrections, 460-640 K). 
For higher temperatures, up to about 720 K (incl. cor- 
rections, 830 K), on-surface adsorbed oxygen is the only 
stable species. For temperatures beyond this, there are 



no stable species except for O atoms adsorbed at under^ 
coordinated surface Ag atoms, such as next to a vacancy 
or step edge, which are stable to 865 K (incl. corrections, 
975 K). From the above analysis, our results point to oxy- 
gen atoms of the thin (4 x 4) surface-oxide structure as 
the main species actuating ethylene epoxidation as this 
occurs experimentally at po 2 = 1 atm. and T = 500- 
600 K, while either low coverages of on-surface oxygen or 
oxygen adsorbed at defects may play a role in the higher 
temperature catalytic reactions. 

In Fig. 3, we consider just the lowest energy struc- 
tures as a function of temperature (T) and pressure (p), 
yeilding the phase-diagram. Here the p-T range in which 
the various structures are stable can quickly be seen. We 
note that the phase-diagram is in-line with available well- 
established experimental results, namely, (i) that under 
UHV conditions, only very low coverage (0.03 ML) of 
chemisorbed O is stable up to about 490 K and (ii) that 
the (4 x 4) phase is stable at low temperatures (and UHV 
pressure) and that it decomposes at 580 Ki2£ Figure 3 
furthermore shows that regardless of the pressure, the 
ordering of the phases is the same. Thus, the same phase 
should be observable at different pressures if the temper- 
ature is adjusted accordingly, which indicates a way to 
bridge the pressure gap - providing that formation of the 
structure is not prevented by kinetics. 

It is interesting to consider the "removal" energy of a 
single (uppermost) O atom of the (4 x 4) structure, which 
is the energy required to move the O atom into the vac- 
uum. This energy is the strongest that we found for 
the O/Ag system, namely, ~3.9 eV, which is due to the 
presence of the lower lying O atoms which stabilizes the 
upper O atoms. Compared to other O/metal systems, 
however, this does not constitute a strong bond. In fact, 
it represents an intermediate value, very similar to the 
identified "optimum" values for highest catalytic activity 
of certain reactions^ Interestingly, recent studies of the 
CO oxidation reaction over Ru have surprisingly found 
that oxides play an active role for reaction, in contrast 
to the hitherto believed pure metali2£i2£ Ruthenium is a 
poor oxidation catalyst since it binds O too strongly, but 
at the high temperature and pressure conditions where 
Ru-oxides form, a new weaker O species appears and the 
system is then catalytically very active. Thus it appears 
that the interaction of oxygen with transition metal cat- 
alysts can serve to "tune" the O-metal bondstrength for 
optimum reactivity through formation of various oxide- 
like species. 

In conclusion, through DFT calculations, and in- 
cluding the effect of the environment, we obtained 
the pressure-temperature phase-diagram for 0/Ag(lll). 
The results show that a proposed molecular ozone-like 
species adsorbed at a surface vacancy is energetically un- 
favorable, as is bulk-dissolved oxygen, which puts a ques- 
tion mark over the latter's hitherto thought important 
role in oxidation reactions. Our results also show that 
silver oxide and oxide-like films can also be ruled out as 
playing an important role due to their low thermal stabil- 
ity. On the other hand, the calculations reveal that silver 
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FIG. 3: Calculated (p, T) phase-diagram for the oxygen- 
Ag(lll) system showing the stable structures. 

can support a thin surface oxide-like structure under the 
temperature and pressure conditions of ethylene epoxi- 
dation and we propose that such atomic O species are 
actuating the reaction. For even higher temperatures, 
only very low concentrations of chemisorbed O are sta- 



ble, or O atoms adsorbed at under-coordinated surface" 
Ag atoms. Therefore these species could possibly play a 
role in the catalytic oxidation reactions that take place 
at high temperature. We note that an understanding of 
a full catalytic cycle requires a kinetic modeling that in- 
cludes all reactant species and intermediates, which, for 
complex systems such as those discussed here, is not yet 
possible. The present investigation represents a crucial 
first step towards this aim, which, using the approach of 
atomistic ab initio thermodynamics, identifies (and ex- 
cludes) oxygen species that are present under the con- 
ditions of catalysis, and helps shed new light upon the 
function of silver as an oxidation catalyst. 

More generally, our results highlight the formation and 
importance of various oxides under catalytic conditions, 
be it in weaker bonded O/noble- metal systems, such 
as O/Ag, or in stronger bonded O/transition- metal 
systems, as e.g., O/Ru. For transition metals that form 
bulk oxides with intermediate heats of formation, e.g. 
PdO, under conditions typical of catalysis, we may 
expect that both bulk oxides and thiner surface-oxides 
could occur, depending sensitively on the pressure and 
temperature. 

We thank D. King and R. Schlogl for helpful discus- 
sions. 
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